Formation of hot heavy nuclei in supernova explosions 
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We point out that during the supernova II type explosion the thermodynamical condition of 
stellar matter between the protoneutron star and the shock front corresponds to the nuclear liquid- 
gas phase coexistence region, which can be investigated in nuclear multifragmentation reactions. 
We have demonstrated, that neutron-rich hot heavy nuclei can be produced in this region. The 
production of these nuclei may influence dynamics of the explosion and contribute to the synthesis 
of heavy elements. 

PACS numbers; 25.70.Pq, 26.50.+X, 26.30.+k, 21.65.-|-f 



In recent years significant progress has been made by 
nuclear community in understanding properties of highly 
excited nuclear systems. Such systems are routinely pro- 
duced now in nuclear reactions induced by hadrons and 
heavy ions of various energies. Under certain conditions, 
which are well studied experimentally, the intermediate 
nuclear system is produced in a state close to statisti- 
cal equilibrium. At low excitation energies this is noth- 
ing but a well known compound nucleus. At excitation 
energies exceeding 3 MeV per nucleon the intermediate 
system first expands and then splits into an ensemble 
of hot nuclear fragments (multifragmentation). At ex- 
citation energies above 10 MeV per nucleon the equili- 
brated system is composed of nucleons and lightest clus- 
ters (vaporisation). These different intermediate states 
can be understood as a manifestation of the liquid-gas 
type phase transition in finite nuclear systems 0. A 
very good description of such systems is obtained within 
the framework of Statistical Multifragmentation Model 
(SMM) 0. 

According to present understanding, based on numer- 
ous theoretical and experimental studies of multifrag- 
mentation reactions, prior to the break-up a transient 
state of nuclear matter is formed, where hot nuclear frag- 
ments exist in equilibrium with free nucleons. This state 
is characterized by a certain temperature T ~ 3 — 6 MeV 
and a density which is typically 3-5 times smaller than 
the nuclear saturation density, po « 0.15 fm~^. Theo- 
retical calculations show that relatively heavy frag- 
ments may survive in the liquid-gas coexistence region. 
In thermodynamic limit these heavy fragments would 
correspond to the infinite liquid phase 0. This statis- 
tical picture of multifragmentation is confirmed by nu- 
merous experimental observations, such as an evolution 
of the fragment mass distribution with energy (tempera- 
ture), fragment correlations revealing the critical behav- 
ior, anomaly in the caloric curve, and by many others 



(see e.g. S IE 0)- Recent experiments (e.g. 'tII) 
directly confirm that primary fragments are really hot, 
their internal excitation energy may reach up to 3 MeV 
per nucleon. Properties of these hot nuclei can be ex- 
tracted from multifragmentation reactions and used in 
analyzes of other physical processes, where similar nuclei 
expected to be produced. 

In this work we are going to use the knowledge accu- 
mulated from multifragmentation studies for better un- 
derstanding nuclear physics associated with collapse of 
massive stars and supernova type II explosions. More 
specifically, we consider the possibility of producing hot 
heavy nuclei in a protoneutron star, and in hot bubble 
between the protoneutron star and the shock front 
This region is crucial for success (or failure) of the super- 
nova explosion. The presence of hot nuclei will influence 
many processes. For example, the electron capture on nu- 
clei plays an important role in supernova dynamics 0. 
In particular, the electron capture rates are sensitive to 
the nuclear composition and details of nuclear structure 
(see e.g. 0). The neutrino- induced reactions are very 
sensitive to the nuclear structure effects and properties 
of weak interactions in nuclei (see e.g. ^3). It is also 
important that the presence of nuclei favors the explosion 
via the energy balance of matter in the bubble [Tjl • 

In the supernova environment, as compared to the nu- 
clear reactions, several new important ingredients should 
be taken into consideration. First, the matter at stellar 
scales must be electrically neutral and therefore electrons 
should be included to balance positive nuclear charge. 
Second, energetic photons present in hot matter may 
change nuclear composition via photonuclear reactions. 
And third, the matter is irradiated by a strong neu- 
trino wind from the protoneutron star. Below we apply 
a grandcanonical version 13] of the SMM to calculate 
mass and charge distributions of nuclear species under 
these conditions. 
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We consider macroscopic volumes of matter consisting 
of various nuclear species {A, Z), nucleons {n — (1, 0) and 
p = (1, 1)), electrons (e~) and positrons (e+) under con- 
dition of electric neutrality. In supernova matter there 
exist several reaction types responsible for the chemical 
composition. At low densities the most important ones 
are: 1) neutron capture and photodisintegration of nuclei 

+ (A+l,Z) + 7, (A,Z) + 7^ (A-l,Z) + 7i, 

, (1) 

2) neutron and proton emission by hot nuclei 

(A, Z) {A~l, T)+n , (A, Z) ^ {A-1, Z-l)+p , (2) 

and 3) weak processes induced by electrons/positrons and 
neutrinos / antineutrinos 

{A, Z)+e- ^ {A, Z-l)+v , {A, Z)+e+ ^ {A, Z+l)+v . 

(3) 

The characteristic reaction times for neutron capture, 
photodisintegration of nuclei and nucleon emission can 
be written as 

Tcap — [{(^nAVnA) Pn] , 
T-^A = [{(^-iAV^a)Pi] , 

respectively. Here ctuA and a-yA are the corresponding 
cross sections, VnA and Vja are the relative (invariant) 
velocities, and r„^p is the neutron (proton) decay width. 
Our estimates show that at temperatures and densities of 
interest these reaction times vary within the range from 
10 to 10^ fm/c, that is indeed very short compared to the 
characteristic hydrodynamic time of a supernova explo- 
sion, about 100 ms [g. The nuclear statistical equihb- 
rium is a reasonable assumption under these conditions. 
However, one should specify what kind of equilibrium is 
expected. For densities p > 10~^po and for the expected 
temperatures of the environment, r< 5 MeV, we obtain 
TyA >> Tcap,Tn,p, i.e. the photodisintegration is less im- 
portant than other processes. There exists a range of 
densities and temperatures, for example, p > lO^^po at 
T = 1 MeV, and p > 10'^ po at T = 3 MeV, where the 
neutron capture dominates, i.e. Tcap < Tn,p. Under these 
conditions new channels for production and decay of nu- 
clei will appear (e.g. a fast break-up with emission of 
a-particles or heavier clusters) which restore the detailed 
balance. We expect that in this situation an ensemble of 
various nuclear species will be generated like in a liquid- 
gas coexistence region, as observed in the multifragmen- 
tation reactions. Here the nuclear system is characterized 
by the temperature T, baryon density ps and electron 
fraction Yg (i.e. the ratio of the electron density to the 
baryon ones). One may expect that new nuclear effects 
come into force in this environment. For example, at high 
temperature the masses and level structure in hot nuclei 
can be different from those observed in cold nuclei (see 
e.g. H). ^ ^ 

The weak interaction reactions are much slower. The 
direct and inverse reactions in eq. Q involve both free 



nucleons and all nuclei present in the matter. It is most 
likely that at early stages of a supernova explosion neutri- 
nos/antineutrinos are trapped inside the neutrinosphere 
around a protoneutron star [T5l |. In this case we in- 
clude the lepton number conservation condition by fix- 
ing the lepton fraction Y^. Out of the surface of the 
neutrinosphere one should take into account the con- 
tinuous neutrino flux propagating through a hot bub- 
ble. Due to large uncertainties in the weak interaction 
rates, below we consider three physically distinctive sit- 
uations: 1) fixed lepton fraction Y^ corresponding to a 
/^-equilibrium with trapped neutrinos inside the neutri- 
nosphere (early stage); 2) fixed electron fraction Y^ but 
without /3-equilibrium inside a hot bubble (early and in- 
termediate times); 3) full /3-equilibrium without neutrino 
(late times, after cooling and neutrino escape). The 
second case corresponds to a non-equilibrium situation 
which may take place in the bubble at early times, be- 
fore the electron capture becomes efficient. 

Chemical potential of a species i with baryon number 
Bi^ charge Qi and lepton number L^, which participates 
in chemical equilibrium, can be found from the general 
expression: 

Pi ^ BifiB + QiPQ + L^Pl (4) 

where /is, fig and pL are three independent chemical 
potentials which are determined from the conservation 
of total baryon number B = J^i Bi electric charge Q — 
Qi and lepton number L = ^ ^ Li of the system. This 
gives 

PAZ = AflB + ZflQ , 

Me- = -Pe+ = -PQ + PL , (5) 
Pi' = -Pi' = PL ■ 

These relations are also valid for nucleons, Pn — PB and 
Pp = Pb + Pq- If and V escape freely from the sys- 
tem, the lepton number conservation is irrelevant and 
fiL = 0. In this case two remaining chemical potentials 
are determined from the conditions of baryon number 
conservation and electro-neutrality: 

PB = TP =Y1 ^P"^^ ' " ^ " X! '^P^^ - Pe^O ■ 

AZ AZ 

Here pe = Pe- — Pe+ is the net electron density. The 
pressure of the relativistic electron-positron gas can be 
written as 

-4(3+ (-)')!. 

Pe \Pe J 

where first order correction due to the finite electron 
mass is included. The net number density Pe and entropy 
density Sg can be obtained now from standard thermo- 
dynamic relations as pe = dPe/dpe and Se = dPe/dT. 
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Neutrinos are taken into account in the same way, but as 
massless particles, and with the spin factor twice smaller 
than the electron one. 

For describing an ensemble of nuclear species under su- 
pernova conditions one can safely use the Grand Canon- 
ical approximation 0, ll.^ . After integrating out transla- 
tional degrees of freedom one can write pressure of nu- 
clear species as 



-Pnuc = ^ X! 9AZ 
AZ 



Vf A^/' 
V A3 



exp 



[FaZ - fJ-Az) 



T^PAz , (6) 



AZ 



where pAz is the density of nuclear species with mass A 
and charge Z. Here gAz is the g.-s. degeneracy factor 

1 /2 

of species {A,Z), At — (27r?i /mArT) is the nucleon 
thermal wavelength, m^r w 939 MeV is the average nu- 
cleon mass. V is the actual volume of the system and 
Vf is so called free volume, which accounts for the finite 
size of nuclear species. We assume that all nuclei have 
normal nuclear density pQ, so that the proper volume of 
a nucleus with mass A is A/ po. At low densities the 
finite-size effect can be included via the excluded volume 
approximation Vf /V « (1 — pb/po)- 

The internal excitations of nuclear species (A, Z) play 
an important role in regulating their abundance. Some- 
times they are included through the population of nu- 
clear levels known for nearly cold nuclei (see e.g. [l6l|). 
However, in the supernova environment not only the ex- 
cited states but also the binding energies of nuclei will 
be strongly affected by the surrounding matter. By this 
reason, we find it more justified to use another approach 
which can easily be generalized to include in-medium 
modifications. Namely, the internal free energy of species 
{A, Z) with A > 4 is parameterized in the spirit of the 
liquid drop model 



Faz [T, Pe 



^ AZ 



Faz 



psym 
^AZ 



^AZ 



(7) 



where the right hand side contains, respectively, the bulk, 
the surface, the symmetry and the Coulomb terms. The 
first three terms are written in the standard form 0, 



^^AziT) = Po 



rp2 _rp2 \ 5/4 



rp2 •• 

— ]A 

£0 



psym 
^AZ 



{A-2Zf 



(8) 

(9) 
(10) 



where wq = 16 MeV, £o = 16 MeV, /3o = 18 MeV, 
Tc = 18 MeV and 7 = 25 MeV are the model param- 
eters which are extracted from nuclear phenomenology 
and provide a good description of multifragmentation 
data 0. S IE la 0- However, these parameters, espe- 
cially 7, can be different in hot neutron-rich nuclei, and 



they need more precise determination in nuclear experi- 
ments (see discussion e.g. in |17|]). In the Coulomb term 
we include the modification due to the screening effect 
of electrons. By using the Wigner-Scitz approximation it 
can be expressed as jl^ 



F' 



AziPe) = gC(pe) 



(11) 



c(Pe) 
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fPe_\ 
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where tq = 1.17 fm and pop — {Z / A)pQ is the proton den- 
sity inside the nuclei. The screening function c{pe) is 1 
at Pe = and at po = Pop- We want to stress that both 
the reduction of the surface energy due to the finite tem- 
perature and the reduction of the Coulomb energy due 
to the finite electron density favor the formation of heavy 
nuclei. Nucleons and light nuclei {A < 4) are considered 
as structureless particles characterized only by mass and 
proper volume. 

As follows from cq. ®, the fate of heavy nuclei de- 
pends sensitively on the relationship between Faz smd 
Paz- In order to avoid an exponentially divergent con- 
tribution to the baryon density, at least in the thermody- 
namic limit (A 00), inequality Faz ^ Paz must hold. 
The equality sign here corresponds to the situation when 
a big nuclear fragment coexists with the gas of smaller 
clusters 3|. When Faz > pAZ only small clusters with 
nearly exponential mass spectrum are present. However, 
there exist thermodynamic conditions corresponding to 
Faz ~ Paz when the mass distribution of nuclear species 
is broadest. The advantage of our approach is that we 
consider all the fragments present in this transition re- 
gion, contrary to the previous calculations [r§.[l9| . which 
consider only one "average" nucleus characterizing the 
liquid phase. 

Mass distributions and charge-to- mass ratios of nuclear 
species were calculated for three sets of physical condi- 
tions expected in protoneutron stars and in supernova 
explosions. We take baryon number B =1000 and per- 
form calculations for all fragments with \< A <1000 and 
0< Z < A in a box of fixed volume V. The fragments 
with larger masses {A >1000) can be produced only at a 
very high density ps > 0.5po, which is appropriate for the 
regions deep inside the protoneutron star, and which is 
not considered here. First we consider the case when lep- 
ton fraction is fixed as expected inside a neutrinosphere. 
Figure 1 shows the results for 1^=0.4 and Yt=0.2. For a 
typical temperature T = 3 MeV we find that the islands 
of heavy nuclei, 100 < A < 500, can appear only at rela- 
tively high baryon density, ps = O.lpo. These nuclei are 
neutron-rich: Z/A « 0.38 and Z/A « 0.27 for Yl=0.4 
and Yj,=0.2, respectively. The Z/A ratios are decreasing 
with A less rapidly than in the nuclear multifragmenta- 
tion case |23|. This can be explained by the screening 
effect of electrons. The width of the charge distribution 
at g iven A is determined by T and 7: tr^ « ^/AT/Wf 
[ij,|23|. At lower density, pB = O.Olpo, the mass dis- 



4 



0,6 



0,5 



0.4 



- 0,3 



0.2 










Yl=0.2 









O.Olp^ 
■ O.OOlpf^ 




250 500 250 

A, fragment mass number 



500 




Y=15MeV 
^-equilibrium 



y=10MeV 

■ O.Olp, 

■ 0.00 Ip^ 

■ O.OOOOlpg 



250 500 250 500 

A, fragment mass number 



FIG. 1: Mass distributions (bottom panels) and average 
charge to mass number ratios (top panels) of hot primary 
fragments produced at temperature T =3 MeV, and for dif- 
ferent densities (in shares of the normal nuclear density po) 
shown on the figure. The calculations are for the lepton num- 
ber conservation: the lepton fraction is Yl =0.4 on the left 
panels, and Yl —0.2 on the right ones. 



FIG. 2: Mass distributions of hot primary fragments pro- 
duced at temperature T =1 MeV, and for different densities. 
Left panels are for the case of the fixed electron fractions Ye 
(no /3-equilibrium) ; right panels are for the full /3-equilibrium 
(without neutrino capture), and for different symmetry en- 
ergy coefficients 7. 



tribution is rather flat up to A « 80 and then decreases 
rapidly for larger A. For p = 10^'^ po only light clus- 
ters are present and the mass distribution drops expo- 
nentially. We have also found that for T = 5 MeV the 
island of heavy nuclei (300 < ^ < 700) is observed at 
very high densities, pB ~ O.Spo, but at pB = O.lpo the 
mass distribution is very broad, up to A w 180. This pic- 
ture is similar to the nuclear liquid-gas coexistence region 
observed in finite systems 

Let us consider the situation more appropriate for a 
hot bubble at early times of a supernova explosion, when 
the neutrino wind from the core interacts with the in- 
falling matter. In this case only the electron fraction is 
fixed, and the electron and proton chemical potentials are 
determined independently, without using the equilibrium 
relation /Ze = —HQ- Corresponding results for Yg = 0.4 
and Ye = 0.2 at T = 1 MeV and several baryon densities 
are presented in Fig. 2 (left top and bottom panels) . One 
can see that heavy nuclei, 50 < A < 500, can be produced 
in very broad range of densities, O.lpo > Pb > lO^^po- 
At given density the mass distribution has a Gaussian 
shape. In the Y^ — 0.4 case the most probable nu- 
clei, corresponding to the maxims of distributions, have 
Z/A ratios 0.400, 0.406, and 0.439, for densities O.lpo, 
10"'^ PO: and 10~^po, respectively. The Gaussian mass 



distributions may in some cases justify earlier calcula- 
tions 0, 01 , when only one heavy nuclear species was 
assumed at each density. As seen from the bottom panel, 
changing the electron fraction from 0.4 to 0.2 leads to a 
significant increase of nuclear masses. Also, the nuclei 
become more neutron rich: the corresponding Z/A ra- 
tios are 0.280, 0.359, and 0.420. It is interesting to note 
that the combination of chemical potentials pe -f /Zp — ^„ 
which drives the electron capture might be quite large in 
this case. For instance, for Ye — 0.4 this value is 14 MeV 
at pb = 10"'^ po, so that the electron capture rate will be 
significantly enhanced compared to the equilibrium rate 
at r = 1 MeVlll. 

In the right panels of Fig. 2 we show results for the 
case of full /3-equilibrium. We observe that the distri- 
butions are quite similar to the ones shown in the left 
panels. For the case of the standard symmetry energy 
with 7=25 MeV we obtain that the Z/A ratios of the 
most probable nuclei are 0.311, 0.359, and 0.435, for den- 
sities O.Olpo, W~^po, and 10~^po, respectively. It is seen 
from the bottom panel that by slightly decreasing the 
7-coe0icient in the symmetry energy one can shift mass 
distributions to higher masses. The nuclei in this case 
are even more neutron rich, the corresponding Z/A ra- 
tios are 0.241, 0.317, and 0.417. For easier comparison, 
in left and right panels we present calculations for the 
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FIG. 3: Average shares of electrons (top panels) and free 
neutrons Yn (bottom panels) versus density. Left panels are 
for the /3-equilibrium with T =1 MeV and different 7. Right 
panels are for the leptoii conservation with T =3 MeV and 
different Yl- 



baryon densities 10~^ and 10~^po, which are more re- 
alistic in the bubble. The comparison shows that in the 
course of /3-equilibration the mass distributions may shift 
by 10-40 mass units. This shift also characterizes uncer- 
tainties in the nuclear composition associated with the 
electron capture reactions. 

Figure 3 shows the fractions of free electrons and neu- 
trons as a function of baryon density. On the left pan- 
els it was calculated for the /3-equilibrium neutrinoless 
matter at T = 1 MeV. These particles, besides heavy 
nuclei, dominate in the supernova matter. As was re- 
alized long ago 0, the electrons are absorbed at large 
densities because of the high electron chemical potential. 
In this work we point out the importance of free neu- 
trons in the liquid-gas coexistence region for maintaining 
a high rate of nuclear reactions. A noticeable change 
in the trend is seen at « 10^"* po- At lower density 
Yn ~ 0.2 that means that 80% of neutrons are trapped in 
nuclei. At higher densities more and more neutrons are 
dripping out of nuclei, and at ps > 10^^ po more than 
half of the neutrons are free. This behavior correlates 
with the decreasing number of electrons and a relatively 
small share of heavy nuclei in the system. However, at 
higher densities, the structure of matter changes because 
of the neutrino/ antineutrino and electron/positron cap- 
ture reactions 12j. The case of lepton conservation is 
shown on the right panels of Fig. 3. One can see that in 



FIG. 4: The mass distribution of heavy elements (with 
A >70) in Solar Sy stem. The experimental data (triangles) 
are taken from [21| [. The lines are calculations at densities 
shown on the figure. The thick solid line is for the lepton con- 
servation at Yl —OA and T =5 MeV. The thin solid, dashed, 
and dotted lines are for the full /3-equilibrium, with the stan- 
dard symmetry energy coefficients 7, and T —1 MeV; the 
dot-dashed line is the same but for the reduced 7. 



this case the number of free neutrons drops with density 
reflecting formation of very big nuclei and transition to 
the liquid phase at ps Po- 

Finally, we find thermodynamic conditions, which al- 
low to produce heavy elements known as the solar cl- 
ement abundances. Here we are not pretending to fit 
the fine structure of these yields, in particular, the pro- 
nounced peaks caused by the structure effects of cold nu- 
clei. Our goal is to propose an explanation for the gross 
distribution of the elements. The results of our fit are pre- 
sented in Fig. 4. It is interesting that by choosing T = 5 
MeV, pb = O.lpo and Yl — 0.4 we are able to achieve a 
reasonable overall fit. This means that elements heavier 
than Fe could be produced in protoneutron stars. How- 
ever, it is difficult to find a mechanism for ejecting this 
material from the star. Another possibility is to consider 
lower densities and temperatures expected in the bubble 
region. In Fig. 4 we also show the fit with T = 1 MeV 
and several baryon densities from 0.003 to lO^'^po- In this 
illustrative calculations we assume that the total baryon 
number is evenly distributed between the regions with 
different densities. A better choice would correspond to 
a temperature and density profile obtained from hydro- 
dynamical simulations of the supernova explosion. As ex- 



6 



pected [Sj, during the explosion, dynamical instabilities 
(convection and other processes) may lead to a large-scale 
mixing of matter and to large density fluctuations. Ap- 
pearance of regions with higher density favors production 
of heavy elements. One can see that we can, in princi- 
pal, explain the element abundances in this way. Also, a 
variation of the parameters of hot nuclei, in particular, 
the 7-coefRcient in the symmetry energy, can influence 
the final predictions. 

One should bear in mind that the mass distribu- 
tions which are presented here correspond to hot pri- 
mary nuclei. After ejection these nuclei will undergo 
de-excitation. At typical temperatures considered here 
(T < 3 MeV) the internal excitation energies are rela- 
tively low, less than 1.0 MeV/nucleon. As well known 
from calculations 2] and nuclear experiments 0, S j 
de-excitation of nuclei with A < 200 will go mainly by 
means of the nuclcon emission. Then the resulting dis- 
tributions of cold nuclei are not very different from the 
primary ones, they are shifted to lower masses by sev- 
eral units. One should expect that shell effects (which, 
however, may be modified by surrounding electrons) will 
play an important role at the de-excitation stage leading 
to the fine structure of the mass distribution. We believe 
that after the de-excitation of hot nuclei, corresponding 
to the time when the ejected matter reaches very low 



densities, the r-process may be responsible for the final 
redistribution of the element abundances leading to the 
pronounced peaks around A ?»80, 130 and 200 |2ll |. 

In conclusion, we have used statistical approach to 
study production of hot heavy nuclei during the collapse 
of massive stars and subsequent supernova explosions. 
Mass and charge distributions of such nuclei have been 
calculated under different assumptions regarding temper- 
ature, baryon density, electron and lepton composition of 
the matter. We have demonstrated that this mechanism 
can contribute significantly to the production of heavy el- 
ements in supernova environment and explain gross fea- 
tures of the element abundances. We also expect that 
the production of hot heavy nuclei may influence the ex- 
plosion dynamics through both the energy balance and 
the capture/production of electrons and neutrinos. 
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